Density modulations in an elongated Bose-Einstein condensate released from a 

disordered potential 
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We observe large density modulations in time-of-night images of elongated Bose-Einstein conden- 
sates, initially confined in a harmonic trap and in the presence of weak disorder. The development of 
these modulations during the time-of- flight and their dependence with the disorder are investigated. 
We render an account of this effect using numerical and analytical calculations. We conclude that 
the observed large density modulations originate from the weak initial density modulations induced 
by the disorder, and not from initial phase fluctuations (thermal or quantum). 
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Gaseous Bose-Einstein condensates (BEC) in disor- 
dered potentials p], 0, H, 0, [B[ offer controllable systems 
to study open basic questions on the effects of disorder 
in quantum media [6]. In this respect, a still debated 
question relies on the nature of the disorder-induced 
superfluid-insulator transitions 0] which can originate 
from strong fluctuations of either the density or the 
phase. This question can be addressed experimentally 
with gaseous BECs in optical disorder p], [H, IIIJH] as 
both density and phase can be measured directly [8[. In 
addition to its fundamental interest, this study is impor- 
tant for BECs on chips [9], and can also shed some light 
on the physics of dirty superconductors [10] and granular 
metals [111 ] . 

In the presence of repulsive interactions, Anderson lo- 
calization is suppressed in a stationary BEC, and weak 
disorder results in small density modulations [12[. How- 
ever, one may wonder whether weak disorder may sig- 
nificantly affect the coherence of a connected BEC and 
entail large phase fluctuations. It has been suggested 
p],[l3[, by analogy with elongated (non-disordered) quasi- 
BECs, that the observation of large fringes in time-of- 
flight (TOF) images of disordered BECs (see Fig. [1] and 
H, [lii [ill ) may signal strong initial phase fluctuations 
Iflllalll. For quasi-BECs, such fringes in TOF im- 
ages are indeed a signature of initial phase fluctuations 
[17l | . However, for disordered BECs, no systematic study 
of these fringes has been reported so far and their relation 
to disorder-induced phase fluctuations is still unclear. 

In this paper, we report a detailed study of the density 
modulations in the TOF images of an elongated, non- 
fragmented three-dimensional (3D) BEC initially placed 
in a weak ID disordered potential. Our main experi- 
mental result is that the fringes in the TOF images are 
reproduced from shot to shot when using the same re- 
alization of the disorder. This excludes disorder-induced 
phase fluctuations (thermal or quantum) in the trapped 
BEC as the origin of the fringes observed after TOF, 
for the parameter range of the experiment (relevant also 
for [l|, ll4|, LL5[ ) . Using analytical and numerical calcula- 
tions, which do not include initial phase fluctuations, we 
show that the fringes actually develop during the TOF 



according to the following scenario. Just after release, 
the initial weak density modulations (induced by the dis- 
order onto the trapped BEC) imprint a phase with axial 
modulations and transversal invariance. Then, the re- 
sulting axial phase modulations are converted into large 
axial density modulations. 

The experiment is detailed in Refs. [2, |5[. We form a 
cigar-shaped BEC of 87 Rb atoms in a Ioffe-Pritchard trap 
of frequencies w z / '2tt = 6. 7Hz and uj_\_/2tt = 660Hz. The 
BEC atom number is No ~ 3 x 10 5 , the length 2Ltf — 
300/im, and the chemical potential ~ 4.5kHz. We 

create a ID speckle (disordered) potential along the z 
axis. The correlation function is C(z) = V^sinc 2 (z/or) 
where both amplitude Vr and correlation length ctr down 
to 0.33/im can be controlled @, [Bj]. The results presented 
in this work correspond to ctr ~ 1.7/im [19]. In the exper- 
iment, we wait 300ms for the BEC to reach equilibrium 
in the presence of disorder and then switch off abruptly 
both magnetic and speckle potentials. We then take ab- 
sorption images of the expanding cloud after a variable 
time of flight £tof w ^ n typical images shown in Fig. [TJ 
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FIG. 1: Upper panel: TOF images of an expanding disor- 
dered BEC for three different times of flight. The vertical axis 
represents the column density along the y axis. Lower panel: 
axial ID density profiles nm(z) (column density integrated 
along the x axis; blue) and ID TF parabolic profiles n? D (z) 
(red). The amplitude of the disorder is 7 = Vr//j, = 0.41. 
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These images show large density modulations along the 
axis z of the disorder. To measure their amplitude, we 
first extract a ID axial density rim(z) by integrating the 
column density over the second transverse direction x. 
We then define rj(z) as the normalized deviations of the 
ID density from the ID parabolic Thomas-Fermi (TF) 
profile n^^z) which fits best the data (red line in Fig. [I]), 
so that uib(z) = n^ B (z)[l + v( z )]- Finally we calculate 
the standard deviation of rj(z) over a given length L: 

Arj = J L dz r] 2 (z) (here, J L dz rj(z) = 0). The calcu- 
lation of Arj is restricted to 70% of the BEC total length 
(L = 1.4Ltf) to avoid the edges where thermal atoms 
are present. Two imperfections reduce the measured At? 
compared to the real value. First, our imaging system 
has a finite resolution L res = 8.5/xm, larger than the vari- 
ation scale ctr of the disorder. This effect is quantified 
by measuring the modulation transfer function (MTF) 
of the imaging system. Second, a slight mis-alignment of 
the probe beam - which is not exactly perpendicular to 
axis z - also reduces the contrast of the fringes [17[ . This 
effect is more difficult to quantify as angles smaller than 
our uncertainty on the probe angle (1°) can drastically 
reduce the contrast. We find that numerics reproduce our 
experimental results assuming a mis-alignment of 0.33°. 

We first study the amplitude of the normalized den- 
sity modulations, At/, as a function of the amplitude of 
the disorder 7 = Vr//i at a given time-of-flight £tof = 
15.3ms (cu±tTOF — 62.2) with experimental results plot- 
ted in Fig. [2^. In the absence of disorder, we observe non- 
vanishing density modulations (At/o — 0.037) larger than 
the noise in the background of the images (Arj n ~ 0.015). 
They are interpreted as small but non-zero phase fluctu- 
ations initially present in our elongated BEC [l6[. The 
calculation of their contribution to the density modula- 
tions in the TOF images as determined in Ref. [17[ agrees 
with our data (see Fig. [2^ and Fig. [3J. In the presence 
of disorder, we find that, for small values of 7 (typically 
7 < 0.2), Arj grows with 7 as Ar) = At/o + 0.64(3)7- For 
larger values of 7, the disordered BEC is fragmented ei- 
ther in the trap or during the expansion, and At? has a 
maximum value of Ar] ~ 0.17 in the experiment [201 ] . 

We also perform numerical integrations of the 3D 
Gross-Pitaevskii equation (GPE) for the expanding dis- 
ordered BEC and extract Arj as in the experiments. The 
numerics do not include any initial phase fluctuations. 
We find a linear dependence of At/ versus 7, At/ c± 3.57 
for the bare numerical results, and Arj ~ 1.237 if we take 
into account the finite resolution of the imaging system 
but not the mis-alignment of the probe (see Fig. (2}d). In 
fact, we find that the numerics agree with the experi- 
ments if, in addition to an offset A?7o to mimic the small 
initial phase fluctuations, we include the systematic cor- 
rection associated to a probe angle of 0.33° (see Fig.[2^). 

We now examine the TOF dynamics of the disordered 
BEC and plot Arj versus £tof m Fig. [3j In the pres- 
ence of disorder, the observed density modulations (red 
points) are clearly enhanced compared to those in the 
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FIG. 2: a) Experimental results (red points) for the density 
modulations at cj^£tof = 62.2 versus the amplitude of the 
disorder. The shaded area corresponds to phase fluctuations 
in our initial elongated BEC, as calculated in Ref. [lTj (the 
error bar reflects the uncertainty on the temperature), b) 
Corresponding numerical results, taking into account the fi- 
nite resolution of our optics but not the misalignment of the 
probe beam. The open blue circles in a) show the same data 
including an offset accounting for the small initial phase fluc- 
tuations and the correction corresponding to the probe angle. 



absence of disorder (blue points). We also observe that 
the density modulations first develop, and then saturate. 
The dynamics of their development is reproduced well by 
our numerical calculations (solid red line), if we take into 
account all imperfections of our imaging system. 

After correcting all experimental imperfections, the 
density modulations in the TOF images turn out to be 
larger (Arj ~ 3.57) than the ones of the trapped BEC 
before TOF (Arj = 27) [2l|. One may wonder whether 
these large density modulations in the TOF images re- 
veal phase fluctuations induced by the disorder in the 
initial BEC p], [l3[ . Actually, several arguments lead us 
to conclude that it is not so. First, our numerics, which 
reproduce the experimental data, do not include initial 
phase fluctuations. Second, the numerical diagonaliza- 
tion of the Bogolyubov equations indicates that the dis- 
order hardly affect the excitation spectrum of a ID BEC 
for the experimental parameters [22[. Last but not least, 
we have observed identical density modulations in suc- 
cessive experiments performed with the same realization 
of the disorder (see also Refs. [2 EH])- Hence, averaging 
over various images taken with the same realization of 
the disorder does not wash out the fringes [23] and this 
excludes initial random fluctuations, quantum or thermal 

We now develop an analytical model for the evolution 
of the BEC density profile during the TOF, which shows 
explicitly how a weak disorder leads to large density mod- 
ulations after a long-enough TOF, without initial phase 
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FIG. 3: Time evolution of the measured density modulations 
A?7 during a TOF for 7 = (no disorder; blue points) and 
7 = 0.41 (with disorder, red points). The shaded area cor- 
responds to the calculation of Ref. [I?} taking into account 
the uncertainty on the temperature. The solid red line is the 
result of numerical calculations for 7 = 0.4 (see text). 



We are thus left with the equations 

dtSh = -(h/m)h d 2 z (2) 
-hd t = g5h/b\ + {h 2 /2m) [\d z 6\ 2 - d 2 5h/2h ]. (3) 

In a first stage, the initial small inhomogeneities of the 
density induced by the ID disordered potential before 
TOF, Sh ~ —V(z)/g [12], hardly evolve since d t Sh(t = 
0) = 0. They are however crucial as they act as an inho- 
mogeneous potential which induces the development of 
a phase modulation 6(z,t) at the beginning of the TOF 
[25(|. From Eqs. (O,©, we find 

0(z, t) ~ arctan(cj ± t) [V(z) /hco^} (4) 
and Sh ~ -V(z)/g - {h d 2 V(z)/mu;l)F(u; ± t) (5) 

where F(r) = j Q T dr f arctan(r / ) =r arctan(r)— In y/l + r 2 . 
From Eq. ([5]), we then find 



fluctuations. Although the probability of fragmentation 
is small for weak-enough disorder, it may happen that 
the BEC is fragmented into a small number of fragments. 
However, for the considered £tof % 1/^, the fragments 
will only weakly overlap as the axial expansion is small. 
Therefore, we neglect fragmentation in our model. 

In the absence of disorder, the TOF expansion of a 
BEC initially trapped in a harmonic potential in the 
Thomas- Fermi regime is self-similar [24[, so that 



A v (t) 



m*)= [11/iW] 1/ %({x j /b j (t)},t)e 



(i) 



with j = 1,..,3 the spatial directions, 0o(r,t) — 
(m/2h) ^Zj(bj/bj)x 2 the dynamical phase and <ft the 
(time-independent) wavefunction of the BEC in the trap. 
The scaling factors bj(t) are governed by the equations 
bj = uo 2 /{bj Y\ k bk) with the initial conditions bj(0) = 1 
and bj(0) = [lij]. In the presence of disorder, we use 
the scaling (pQ) and we write the (now time-dependent) 
wavefunction <j)(p,z,t) = y/h(p, z, t)e l9 ^ p,z,t ^ where p = 
^x 2j ry 2 is the radial coordinate. In the absence of phase 
fluctuations, </>(p, z, t = 0) is real (up to a homogeneous 
phase) as it is the ground state of the trapped, disordered 
BEC. The TOF dynamics is then governed by two cou- 
pled equations for the density ft and the phase 9 which 
are equivalent to the complete time-dependent GPE. 

Let us introduce now a couple of approximations. 
First, in elongated 3D BECs, the expansion for uj z t < 1 is 
mainly radial and b z (t) ~ 1. Second, we assume not too 
large perturbation of the density, so that h = ho+5h with 
no, the density in the absence of disorder and Sh <C ho. 
Using local density approximation, we can neglect all 
the spatial derivatives of ho. We also neglect the radial 
derivatives of 5h since the ID disorder induces short- 
range spatial inhomogeneities mainly along the z-axis. 



p 



2 

3 V hw± 



2 / ^ x 2 



. (6) 



Hence, A77 first slightly decreases. It is easily understood 
as an interacting BEC initially at rest will tend to fill its 
holes when released from the disordered potential. The 
solution dU),© is valid as long as the contribution of the 
last two terms in Eq. (j3]) remains small [i.e. for co±t <C 
(cr/O 2 an d u±t <C (ctr H){fujj±_ j V^raO] • I n addition, it 
requires that the density modulations do not vary much 
(i.e. the second rhs term in Eq. ([5]) is small compared to 
the first one). For the experimental parameters, the last 
condition is the most restrictive. It defines a typical time 
t = (l/co±)F- 1 [{o- R /£) 2 (hu; ± /p) 2 /2} during which the 
radial expansion imprints a phase modulation due to the 
initial inhomogeneities of the BEC density created by the 
disorder before the TOF. In particular, if to ^> the 
phase modulations freeze at 6{z) (jr/2)[V(z)/huj±\. 

In a second stage, the phase modulations are con- 
verted into density modulations similarly as thermal 
phase fluctuations do during the TOF of an elongated 
quasi-BEC [l7| • For t > t\ where t\ is a typical 
time much longer than l/u± (see below), the scal- 
ing parameter b±(t) becomes large so that the first 
rhs term in Eq. ([3]) can now be neglected. Assum- 
ing small phase gradients, we are left with the equa- 
tion d 2 Shk + h 2 k A 5hk/ '4m 2 = where 5hk(p,t) is the 
ID Fourier transform of 5h along z, and whose so- 
lution reads Shk(t) = Shk(ti) cos [(hk 2 /2m)(t — ti)] + 
(2m5h k (t 1 )/hk 2 )sm [(hk 2 /2m)(t - h)]. If t > 
we can take 1/uj± -C t\ < to and the exact value of t\ 
does not matter much (we use t\ = to). If to ^ l/^±, 
the determination of t\ is not straightforward but can 
be found through fitting procedures, for instance. Then, 
according to Eq. ([5]), b~hk(t\) c± —V(z)/g is mainly deter- 
mined by the initial density modulations of the trapped 
BEC while 5hk(ti) — — (hod 2 V(z)/mu±) arctan(o;^ti) 
results from the phase modulations created in the first 
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FIG. 4: Dynamics of density modulations as obtained numerically (solid red line with error bars) and comparisons with 
Eq. J6} (green dashed line) and Eq. (0 (blue dotted line), a) The parameters are the same as in the experiment (in particular, 
ctr = 1.7/xm) with Vr = 0.02/i. Here u±to ~ 3 and we have used t\ — 0.5£o- b) Same as a) but with ctr = 3.4/xm. Here, 
oo±to c± 8.5 ^> 1, so that we have used t\ = to (see text). 



stage of the TOF. For hou± <C \i as in the experiment, 
the cosine term can be neglected and we find 

Ar](t) ~ Vs(Vr/Huj±) arctan(cj ± ti)/[cr R , t - t{\ (7) 

where J(cr R ,£) = y dtz (1 — k) sin 2 [(2ht / ma^) k 2 } for 
a speckle potential. 

Numerical integrations of the complete 3D GPE con- 
firm the expected behavior of An during the TOF at 
short [Eq. (j6])] and long [Eq. Q] times as shown in Fig.JH 
This validates our scenario in a quantitative manner. 

Three remarks are in order. First, we find that, due to 
the development of phase modulations in the first stage of 
the TOF, the density modulations in the expanded BEC 
(An oc Vr/Hujj_) can be larger than those in the trapped 
BEC (A77 oc Vr///) if ji > hco±. Second, the density pat- 
tern is completely determined by the realization of the 
disorder. Third, Eq. (0) shows that the density modu- 
lations saturate at An ~ V^(Vr/Huj±) arctan(o;_Lti) for 
very long times t [26] . These properties are in qualitative 
agreement with the experimental observations. 

In conclusion, we have shown that the large fringes ob- 
served in this work and Refs. [2, [3, [l5[ m TOF images of 



disordered BECs do not rely on initial disorder-induced 
phase fluctuations. They actually result from the TOF 
process of a BEC with small initial density modulations 
but without phase fluctuations. A phase modulation de- 
termined by the weak initial modulations of the BEC 
density first develops and is later converted into large 
density modulations. Our analytical calculations based 
on this scenario agree with numerical calculations and 
experimental observations. Nevertheless, our results do 
not exclude that, in different regimes, disorder might en- 
hance phase fluctuations. This question is crucial in con- 
nection to the nature of superfluid-insulator transitions 
in the presence of disorder. Dilute BECs can help an- 
swering it as their phase coherence can be probed with 
accuracy [17|, |_18f] . Revealing this possible effect via TOF 
images requires taking into account the phase modula- 
tion which develops during the first stage of the TOF as 
demonstrated here. 
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